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1 Introduction 



The Lifshitz black holes [l]l2i[3j|IJ[5j[6j[7J|8] have received considerable attentions since 
these may provide a model of generalizing AdS/CFT correspondence to non-relativistic con- 
densed matter physics as the Lif/CFT correspondence [UEBEI]- Although their asymptotic 
spacetimes are known to be Lifshitz, it is a non-trivial task to find an analytic solution. 
One of the known solutions is a four-dimensional topological black hole which is asymp- 
totically Lifshitz with the dynamical exponent z = 2 [12]. Analytic black hole solutions 
with planar horizon were found in the Einstein-scalar-massive vector theory [13] and in the 
Einstein-scalar-Maxwell theory [14]. Another analytic solution has been recently found in the 
Lovelock gravity [15]. The z = 3 Lifshitz black hole [16] was derived from the new massive 
gravity (NMG) [17]. Numerical solutions and thermodynamic property of Lifshitz black hole 
were explored in [HI [19] . 

A thermodynamic study is important to understand the Lifshitz black hole because heat 
capacity and free energy determine the global stability of the Lifshitz black hole. A posi- 
tive (negative) heat capacity imply thermally stable (unstable) black hole, while a positive 
(negative) free energy means unfavorable (favorable) configuration in given ensemble. Hence, 
a black hole with positive heat capacity and negative free energy is considered as a globally 
stable black hole ( GSBH). However, the thermodynamic study on Lifshitz black holes was 
limited because it was difficult to compute their conserved quantities in Lifshitz spacetimes. 
Recently, there was a progress on computation of mass and related thermodynamic quanti- 
ties by using the ADT method [20j [21] and the Euclidean action approach [22]. Concerning 
the mass of 3D Lifshitz black hole, there is an apparent discrepancy between J\A — 8G3I 4 
obtained from the ADT method [20] and M. = 4G ^ 4 from other methods [221 [23], [21] . Phase 
transitions between Lifshitz black holes and other configurations were investigated by using 
on-shell and off-shell free energies [25J. 

On the other hand, quasinormal modes (QNMs) of a perturbed field contain important 
information about the black hole. Their quasinormal frequencies (QNFs) are given by u> = 
lor — iui whose real part represents the perturbation oscillation and whose imaginary part 
denotes the rate at which this oscillation is damped, because of the presence of black hole 
horizon. In this sense, one requires ui > which is consistent with the stability condition of 
the black hole. Since in asymptotically AdS (Lifshitz) spacetimes, spacelike infinity acts like 
a reflecting boundary, Dirichlet, Neumann, or mixed boundary condition have to be imposed 
there. The QNFs could be obtained by solving the Klein-Gordon equation for a minimally 
coupled scalar by imposing the boundary conditions: ingoing mode near the horizon and 
Dirichlet condition at infinity. 
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Importantly, all known QNFs ^ZE\ [27J [2S] have no real part, which implies that they are 
purely imaginary. If the purely imaginary frequency (cj = —iuj, coj > 0) represents an 
interesting feature of the Lifshitz black hole, it is very curious to explore its connection to 
thermodynamic property We note that the imaginary part of QNFs involves the temperature 
of black hole and a GSBH may provide an analytic form of QNFs. 

According to the Hod's conjecture [29], the asymptotic QNFs is related to the quantized 
black hole area. Identifying the vibrational frequency lu(E) with the real part ur, it leads to 
an area quantization of AA n = 41n[3]£p which is not universal for all black holes. For a large 
damped case, Maggiore [30] has proposed that the identification of co(E) with the imaginary 
part uj might lead to the Bekenstein universal quantization of AA n = 87r£p [3Tj . Hence, the 
analytic computation of QNFs is crucial for extracting an important information on Lifshitz 
black holes obtained from different gravitational theories. 

In this work, we investigate Lifshitz black holes by exploring the connection between ther- 
modynamic property and quasinormal frequencies. It is shown that globally stable Lifshitz 
black holes provide purely imaginary quasinormal frequencies. We also find that by using 
the Maggiore's method, the horizon area and entropy can be equally spaced for these black 
holes. 

The organization of our work is as follows. In section 2, we study the 3D Lifshitz black 
hole by exploring a connection between its thermodynamic quantities and QNFs, where 
QNFs are already known in [26]. We investigate the 2D Lifshitz black hole by obtaining its 
purely imaginary QNFs in section 3. In section 4, we study two 4D Lifshitz black holes: one 
obtained from the Einstein-scalar-massive vector theory and the other from the Einstein- 
scalar-Mawxwll theory. QNFs of the former black hole were found in f27\ ¥2E\ when replacing 
a radial coordinate r by 1/r, while we obtain newly QNFs of the latter. Finally, we find the 
area spectrum of two 4D Lifshitz black holes after reviewing the 3D Lifshitz black hole in 



The NMG [T7] composed of the Einstein gravity with a cosmological constant A and higher- 
order curvature terms is given by 




2 3D Lifshitz black hole 
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where G3 is a three-dimensional Newton constant and fh 2 a parameter with mass dimension 
2. We mention that to avoid negative mass and entropy, it is necessary to take "— " sign in 
the front of [S^ H + Sjg C ] . The field equation is given by 

- -g^Tl + Ag^ - t^K^ = 0, (4) 



where 



+ m^pa-RT --nn tlv -n pa n pa g flu + -n 2 g flv . (5) 

Z o 

We have to choose fh 2 = — 7^2 and A = — S to obtain the z = 3 Lifshitz black hole solution. 
Here i the curvature radius of Lifshitz spacetimes. Explicitly, we find the Lifshitz black hole 
solution [TB] as 

tZalo = ^d^tfe" = - + ZT]^ + r ^ ' (6) 

where M is an integration constant related to the the mass of black hole. From the condition 
of g rr = 0, the event horizon is determined to be r = r + = tyM. This line element is 
invariant under the anisotropic scaling of 

r M 

t^XH, 0^A0, r->-, M^—. (7) 

r 2 

For z — 1, the ADM mass is determined to be M = while for z = 3, the ADM mass is 
not yet fixed completely. 

All thermodynamic quantities were driven by using the Euclidean action approach. Its 
Hawking temperature (T#), mass (/A), heat capacity (C = J^- ), Bekenstein-Hawking en- 
tropy (Sbh), and Helmholtz free energy (F = M. — T h Sbh) are given by 

T _ r + w _ r + r _ 47rr + o _ 27rr + p _ 3r + /ox 

At this stage, we mention a global structure of Lifshitz black hole. Its Penrose diagram 
is figured out to be M where a light-like curvature singularity is located at r = (top 
and bottom), while Lifshitz asymptote is at r = 00 (two sides). Generally, a black hole is 
increasing by absorbing radiations in the heat reservoir, while a black hole is decreasing by 
Hawking radiation as evaporation process. In studying the phase transition, two important 
quantities are the heat capacity C which shows thermal stability (instability) for C > 0(C < 
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0) and free energy F which indicates the global stability for F < 0. For the case of positive 
heat capacity and negative free energy, we call it the globally stable black hole (GSBH). Here 
it is observed from (jHJ) that the Lifshitz black hole belongs to the GSBH because of C > 
and F < 0. 

In order to make a connection to the QNMs, we consider the minimally coupled scalar 
described by the Klein-Gordon equation 



□3D - m 



(p = 



(9) 



in the background of Lifshitz black hole ([H]) for z = 3. Decomposing ip with y = r + /r as 

v(t,y,<P) = R(y)e- MK *, (10) 
the radial equation takes the form 

y 2 



R" + 



■ 2 3 



2/(1 - y 2 ) 



R' + 



i-y 2 



u 2 y* 



m 
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R = 0. 
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Here we note that r e [r + , oo) is mapped inversely to y G [1, 0). The solution to (ITT]) is given 
by the confluent Heun (HeunC) functions as 



R( y ) = C iy 2+a {l -i/ 2 )2HeunC 
+ C 2 y 2 - a (l -y 2 )^ HeunC 
where C\ and Ci are arbitrary constants and 
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a = 2\/l + —— > 2, /3 



0, -a, (3, - 
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4 ' 4 4M 



:i2i 



2nT, 



(13) 



Requiring the Dirichlet condition at infinity (y = 0) leads to C 2 = because of 2 — a < 0. 
In order to impose the ingoing mode at horizon, one uses the connection formula 



HeunC 



0,a,/3,7,5; z 



T(a + 1)T(-P) 
T(l — (3 + K)T(a — K) 



HeunC 



0,/3,a, -7,7 + 5;l-z 



+ 



r(« + i)r(/3) 



■(1 - ^) _/3 HeunC 0, -/?, a, -7, 7 + 8; 1 - z 



r(l + p + S)T(a - S) 
where K and S are determined by solving two algebraic equations 



(14) 



K 2 + (1 - a - 0)K - a - f3 - e 
S 2 + (1 - a + P)S - a - a(3 - e + 
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with e = [1 - (a + + l)]/2 - 5. Near the horizon (y -» 1), using PU), the solution flTJ) 
can be written by 



Ci 



where £i and £ 2 are given by 

r(i + a )r(-/3) 



*i(i-i/ 2 
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r(i + «)r(/3) 



r(a - it)r(l - /9 + A") ' ^ z r(a-5)r(l + /3 + 5) 

For the real u = lor in f|T5|) . the former (latter) correspond to ingoing mode 
mode <( — \y=i) because the scalar field cp behaves as 



<P 



tie^l - y 

£ie 2wT ' 



+ 6e 



" + 6e 
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Z > 2 

In(l-y)] 



2ttT 



(15) 
(16) 

| y =i (outgoing 
(17) 



near the horizon. Hence, to obtain the ingoing mode at the horizon, one imposes £2 = 
which implies that T(a — S) — > 00 or T(l + (3 + S) — > 00. This could be done by requiring 



a — S = —n, 1 + (3 + S = —n, n 

which lead to the same expression of quasinormal frequency. 
Finally, one recovers the purely imaginary frequency 



0,1,2, 
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l-2n-(4+m 2 1/2 ±(7+ 



3m 



2 02 



K 
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+ {3+6n){4+m 2 e 2 ) 1/2 +6n{n+r 
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Its asymptotic frequency takes the form 



H 



n. 



(19) 



(20) 



when choosing "o; + " because selecting "o;_" leads to unstable quasinormal modes. At this 
stage, we note that the expression ffT9|) has already appeared in [26]. However, it seems that 
in deriving this expression, they have made two mistakes: one was to choose C\ = 0, instead 
of C2 = 0, and the other was to use a wrong connection formula, instead of a correct one 
(JSJ). Two mistakes happen to provide a correct expression ([T9"j) . In comparison to the BTZ 
quasinormal frequency 



.BTZ 
± 



±- 



iAttT, 



H 



n+ -(1 + v 7 ! + m 2 £ 2 ) 



(21) 



the imaginary frequency involves the angular quantum number k in the Lifshitz black hole. 

Consequently, we confirm that the QNFs of 3D Lifshitz black hole, which is a GSBH, are 
purely imaginary. We check the connection between purely imaginary quasinormal frequency 
and positive heat capacity and negative free energy for the 3D Lifshitz black hole. 
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3 2D Lifshitz black hole 



It is necessary to study a lower dimensional Lifshitz black hole because the exact computation 
of QNFs is available for the lower dimensional gravitational system. The only known 2D 
Lifshitz black hole could be obtained when applying the Achucarro-Ortiz (AO) dimensional 
reduction to the NMG action 



ds 2 {3) = g ij dx i dx j + £ 2 $ 2 d6 2 (22) 

with the dilaton <£>. Integrating over 9 on S 1 , the action (JT|) reduces to the 2D effective dilaton 
action as 



Snm 



NMG 



Seh + Shc 



(23) 



where 



S 



EH 



Shc 



J d 2 aV=^$(i? (2) -2A), 



2m 2 



d 2 xy/^g $ 



Shc contains fourth-order derivatives as the dilatonic kinetic term. It turned out that for 
$ = rji and z = 3, equations of motion for 2D metric tensor g l i and dilaton $ admit the 2D 
Lifshitz black hole solution 



^^2D — CJijd$ 



2\ 3 



1 



M£ 2 \ , 2 dr 2 
dt 2 + - 



M 



(24) 



All thermodynamic quantities of 2D Lifshitz black hole are the same as (|S]) of the 3D Lifshitz 
black hole because the AO-dimensional reduction preserves all thermodynamic properties of 
3D Lifshitz black hole. Hence, the 2D Lifshitz black hole is also a GSBH. 
Now we introduce a minimally coupled scalar equation 



□on - rn 



in the background of 2D Lifshitz black hole 

y = i^Mjr = r + /r as 

p) = Me** 



= (25) 
to find the QNMs. Decomposing ip with 

(26) 



the radial equation becomes 



P"- 



p' + 



m 



M 3 (l-y 2 



p = 0, 



(27) 
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which is similar to the 3D radial equation (ITT]) with k = (s-mode). Also we note a coordinate 
mapping: r G [r + , oo) — >■ y G (0, 1]. Solving the equation (|2"T|) leads to the solution which is 
expressed in terms of the HeunC functions as 



p(y) = C W l^(l-y 2 )^ReunC[o n ,{3,-^?l + ^y 2 
+ C 2 y^\\ - y 2 )^HeunC [o, -7, A -y , ^ + ^; y 2 
where C\ and C*2 are arbitrary constants and 



(21 



3 / 4m 2 £ 2 3 „ . u£ 
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— 



27rT, 



(29) 



We observe that at infinity (y — > 0), imposing the Dirichlet condition leads to C 2 = 0. In 
order to obtain the ingoing mode at horizon, we use the connection formula (fl4|) . Near the 
horizon (y — > 1), the solution ( 128|) takes the form 



Py^i — C\ 



6(1-^+6(1-^)^ 



r(i + 7 )r(/3) 



where the coefficients £i and £ 2 are given by 

~ = r(i + 7 )r(-/3) ~ 
?1 r( 7 - iqr(i - + ' r( 7 -s)r(i + /3 + 5) 

In these expressions, K and S are determined by solving two equations 
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(30) 



(31) 



K 2 + (1 - 7 - P)K - 7 - 13 - e + 
5 2 + (l-7 + /3)5- 7 - 7 /3-e + 



5 



0. 







with 5 = -/3 2 /4 and e = [1 - (7 + l)(/3 + l)]/2 - 7 2 /4 - 3/16. 

For the real u; = in (1301) . the former (latter) correspond to ingoing mode (outgoing 
mode) because near the horizon, the scalar field ip takes the form 



= £ l6 — *(i-y^ +^-^(1-^-2 



(32) 



To obtain the ingoing mode at the horizon, one has to impose £2 = which means that 
I\7 — S) — > 00 or T(l + {3 + S) — > 00. This is achieved by requiring 



7 - £ = -n, 1 + + S = -n, n = 0, 1, 2, 



(33) 
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which lead to the same QNFs. 

Consequently, one obtains the purely imaginary frequency 



AttT, 



H 



/9 o f o\ 1/2 /19 3m 2 £ 2 , ,/9 ,,\i/2 . .\i/2 
-l-2n-[-+m 2 £ 2 ) ±[ — + ^— + {3 + 6n)[-+m 2 e 2 ) +6n(n + l)J 

(34) 

We note that its asymptotic frequency takes the same form as that of 3D Lifshitz black hole 
in flUD 

3 -±^ = -i(V6-2)n (35) 

which shows that 2D Lifshitz black hole is very similar to 3D Lifshitz black hole. In compar- 
ison, we introduce the QNFs for the AdS2 black hole 



-i\n-\ (36) 



2ttT, 



H 

which is purely imaginary. 

As a result, it is shown that a scalar propagating in the 2D Lifshitz spacetimes provides the 
purely imaginary quasinormal frequency. We find the connection between purely imaginary 
quasinormal frequency and positive heat capacity and negative free energy for 2D Lifshitz 
black hole. 

4 4D Lifshitz black holes 

Up to now, we have found the connection between purely imaginary quasinormal frequency 
and positive heat capacity and negative free energy for 2D and 3D Lifshitz black holes. In 
order to confirm this suggested connection, it is necessary to investigate 4D Lifshitz black 
holes. Two known Lifshitz black holes were obtained from the Einstein-scalar-massive vector 
theory and the Einstein-scalar-Maxwell theory. 

4.1 Einstein-scalar- massive vector theory 

First we study the 4D1 Lifshitz black hole obtained from the Einstein-scalar-massive vector 
theory [13] 



Sesmv = tttV / d 4 xV=g[R -2A- m 2 A^ - 2(e^ - 1) - \e~^F^F^ (37) 

107TG4 J L 



where 



A = - - +z + 4 ; m 2 = 2z, F — dA (38) 
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with L 2 = 1. The Lifshitz black hole solution for z = 2 is given by 

dr 2 



U6 4D1 



-20 



-r 2z f(r)dt 2 + r 2 (dx\ + dx\)j + 
f(r)r 2 dt 



f(r)r 2 



1 + A 



with the metric function 



/(r) 



V2 



„2 ■ 



(39) 



It is important to note that the above solution could be obtained by replacing r by 1/r in 
the original solution appeared in pQ. 

Applying the Euclidean action approach to this theory, one finds thermodynamic quan- 
tities 

- 2 V 2 ri „ V 2 rl n r\V 2 „ V 2 r A , 



T * = £> M = 16^? C = 4G? ^ 4G? F = -T^> ^ 
where V 2 is the volume of the transverse directions. The 4D1 Lifshitz black hole is a GSBH 
because of C > and F < 0. We note that the thermodynamic quantities in f l40i) are 
obtained by replacing r + by l/r+ in the original expressions as 

1 v 2 „ v 2 v 2 v 2 



H 



M 



C 



s 



lTir 2 + 167rG 4 r£' " 4G 4 r 2 ' " BH 4G 4 r^' x 167rG 4 r£" 

Here, the limit r + — >■ of all thermodynamic quantities goes to infinity which indicates that 
the original radial coordinate is not appropriate for describing thermodynamics of the Lifshitz 
black hole. In addition, the computation of QNFs has been performed by using the original 
coordinate appeared in [27], showing four types of QNFs: Ui, i — 1, 2, 3, 4. This is the reason 
why we use the new radial coordinate as the inverse of original coordinate. 

In order to compute QNFs for the z = 2 Lifshitz black hole (I3"9~j) . we consider a massive 
scalar field given by 



(41) 



□ 



4D1 



m 



tfi = 0. 



Assuming = H\(r)e %wt e l ( k ^ x i+ k ^) anc j introducing y = r + /r, equation 



y- 



-H'i 



yW 



r Ar\m + k y )(1 — y ) 



,2\2 



#i = 0, 



(42) 
becomes 

(43) 



i/(i - v 2 Y~ L ' ^+y 2 (i - y* 

where k 2 = kf + k\ and the prime (') denotes the differentiation with respect to It is found 
that the solution to Eq.f l43|) is given by the hypergeometric functions 



EM = D iy a (l- y y a F a 



a, 6, c; y 2 



+ D 22/ 4 -«(l - 2/Y 2 F, 



a — c + 1, b — c + 1, 2 — c; 



(44) 
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where 2F1 is the hypergeometric function, Di and D 2 are arbitrary constants, and 

a = 2 



1 + A 1 



n? 



2n 



>4, /3 



1UJ 



AttT, 
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In the hypergeometric function, the arguments of a, b, and c are given by 



a 1 1 



2 2 2rj 



ri — u 2 — r+k 2 , 



a 1 1 
^2-2 + H 
a — 1. 



2 JL2 



(45) 



(46) 

(47) 
(48) 



At infinity (y — > 0), we see that D2 of (jUJ) should be zero because it corresponds to non- 
normalizable mode (4 — a < 0). At the horizon (y = 1), using the connection formula 
which connects y 2 to 1 — y 2 for the hypergeometric function, the solution (1441) becomes 



ffift/) = D iy a 



1 _ ^ r(c)r(c-a-6) 



r(c-a)r(c-6) 
r(c)r(a + 6-c 



r(o)r(6) 



2^1 



>-Fi a, 6, a + b — c + 1; 1 — y 2 
c — a, c — 6, c — a — b + 1; 1 — y 2 ] 



(49) 



It is worth to note that at the horizon, the first term in fT49|) corresponds to the ingoing 
mode, while the second corresponds to the outgoing mode. The ingoing mode near the 
horizon (y — > 1) is obtained by imposing the condition 



T{c)T{a + b-c) 



0. 



which implies that 



-n, b = —n, n = 0, 1, 2, • • • . 



(50) 



(51) 



Consequently, we find the purely imaginary quasinormal frequency as 

'(2n + 1) {{An 2 + An - 4 - m 2 ) + ^) 



^4D1 



27TT; 



—I 



H 



2(4n 2 + 4n-3-m 2 ) 



+ 



((4n 2 + 4n - 2 - m 2 ) - ^4 + 

2(4n 2 + 4n - 3 - m 2 ) 



(52) 
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which is the same expression as those of U2 in [27] and u in [2S]. The n = case leads 
to u> j > 0, which means that the 4D1 Lifshitz black hole are stable against the external 
perturbations. Also, we observe that in the large n limit, CU4D1 becomes 



y 4Dl 



2ttT, 



H 



-in, (53) 



which will be used to derive the area quantization in section 5. 

Finally, we confirm that a scalar propagating in the 4D1 Lifshitz black hole has purely 
imaginary quasinormal frequency. We find the connection between purely imaginary quasi- 
normal frequency and positive heat capacity and negative free energy for the 4D1 Lifshitz 
black hole obtained from the Einstein-scalar-massive vector theory. 

4.2 Einstein-scalar-Maxwell theory 

We introduce the action of Einstein-scalar-Maxwell theory [T4] 

Sesm = Y^Gl j d * x ^9l R - 2A " \ d ^<t> ~ \e^F, u F^l (54) 

where A is the cosmological constant and two fields are a massless scalar and a Maxwell field. 
It admits the 4D2 Lifshitz black hole with z = 2 as solution to equations of motion 



2 



<7s 2 - T 2 

U6 4D2 — ^ 



1=1 



r z + 2 ~ 1 4 
/(r) = l-^, ^ = 1, A 2 ' 



t H -qr , A — — , 

q 2 = 2L 2 (z~l)(z + 2) } (55) 

where the event horizon is located at r = r+. This line element is invariant under the 
anisotropic scaling of t — > X z t,Xi — > \Xi,r — > r/A, and r + — > r + /A. It is important to 
note from the last relation of (155]) that the charge q is not an independent hair because it is 
determined by the curvature radius L of Lifshitz black hole and its dynamical exponent z, 
which contrasts to the Reissner-Nordstrom-AdS black hole. A similar case was found in the 
charged MTZ black hole [3TI 13"%] 

The temperature and Bekenstein-Hawking entropy are determined by 



rpz 



: "V;, S BR = ^rl, (56) 



4tt 
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where V 2 denotes the volume of two-dimensional spatial directions. Mass, heat capacity, and 
Helmholtz free energy are obtained by using Euclidean action approach as 



M z 



2L 2 V 2 r z + 
16vrG 4 



2+2 



C z 



dM z 



2L 2 V 2 r\ 
AzGa ' 



zL 2 V 2 r z + +2 
16ttG 4 '' 



(57) 



16ttG 4 



-. We stress that 



where F z is different from the Gibbs free energy defined by F z = 
the 4D2 Lifshitz black hole is also a GSBH because of C z > and F z < 0. 

In order to compute QNFs for the z = 2 Lifshitz black hole (155]) . we consider a massive 
scalar field given by 



□402 - m 



#2=0. 



(58 



Assuming #2 = H 2 (r)e %u)t e H^izi+faa^) anc [ introducing y = r + /r, equation fl58l) becomes 

(59) 



-Ho - — -n 9 H 3 — ttt; — n 2 = U, 



2/(1 — 2/ 4 



r^y 2 (l - y 



4^2 



where k 2 = k\+k 2 and the prime (') denotes the differentiation with respect to y. Importantly, 
the solution to fl59l) is given by the general Heun (HeunG) function as follows: 



D iy 2+a (l 



H 2 (y) 



y 2 y*(l -jy^HeunG 



ft 



1, -1, 1 + 1 + 1 + a, l + (3; y 2 



2 - a ^ ' ~ 2 ^-f(l -y 2 )^ HeunG 



ft 



7 + 2«/3, 1- -, ] 



a 
2' 



where and D 2 are arbitrary constants, and 



m 2 L 2 

a = 2\ 1 + — — > 2, 



/3 = -i 



00 



2nT, 
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ft, 1 + /3; y 2 
(60) 



(61) 



47rTi 



At infinity (7/ = 0), we see that D 2 should be zero because it corresponds to non-normalizable 
mode (2 — a < 0). In order to derive the asymptotic form of the HeunG function near the 
horizon (y — > 1), we use the formula [391 1^0] 



HeunG 



61, b 2 , ci\, a 2 , CI3, a 4 ; z 



HeunG 



bi, b 2 — (a 4 — 1)0361, a 2 — a 4 + 1, ai — a 4 + 1, a 3 , 2 — a 4 ; z 



E\ HeunG 



1 — 61, — b 2 — a\a 2l ai, 02, 1 + ai + a.2 — 03 — a 4 , a 4 ; 1 — z 



+ 



' 1 - ^)«3+a 4 -ai-a 2jE2 ReunG 



1 — bi, —b 2 — a x a 2 — (a 3 + a 4 — a x — a 2 )(a 3 + a 4 — 6xo 3 ) 



03 + a 4 — Oi, a 3 + a 4 — ^2; 1 — cti — 02 + 03 + a 4 , a 4 ; 1 — z 



(62) 
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where E\ and E 2 are given by 



61, b 2 , Oi, a.2, 04; 1 



61, b 2 — bia 3 (d3 + a 4 — di — CL2) , 03 + 04 — ai, 03 + 04 — 02, a 3 , a 4 ; 1 



(63) 
.(64) 

(65) 



We note that the first term in fl65l) corresponds to the outgoing mode, while the second 
corresponds to the ingoing mode. To obtain the ingoing mode for the solution (J65]) near the 
horizon (y — > 1), we impose E\ = which implies 



Ex = HeunG 
E 2 = HeunG 

Using the formula fl62j) . the solution fl60|) takes the form near the horizon 
ff 2 (y) = D iy 2+a (l + y 2 fA{l- tfY-Ex + (1 - y 2 ) W • 



HeunG 



-1, - 



k 2 



4r 



2 ■ 



1 + - - 0, 1 + - - 0, 1 + a, 1 - 0; 1 



0. 



(66) 



At this stage, we have to mention that a little bit of general Heun function and its connec- 
tion formula is known, in comparison with the HeunC function and hypergeometric function. 
Hence, it seems to be a formidable task to obtain QNFs from a general condition of (166]) . 
Fortunately, we observe that for k 2 = (s-mode), f)66p reduces to a condition for the hyper- 
geometric function 



1 a. 8 1 a 8 a 

- H -, - H -, 1 + -; 1 

2 4 2' 2 4 2' 2' 



In deriving this condition, we used the reduction formula [4Tj 



HeunG 



0-1 + a 2 — a 3 + 1 

1, 0, ai, ci2, 03, ; z 



i F 1 







ax a 2 1 + a 3 2 

T' T' 2 ' 2 . 



(67) 



(6£ 



Considering a relation for the hypergeometric function 

Alcu c, c 3 ; 1] = r(c3)r(c 3 - c 1 -c 2 ) 

T(c 3 - Ci)r(c 3 - C 2 ) 

the condition of 2 E± [ci, C2, C3; 1] = implies that 

C3 — ci = — n, or C3 — C2 = —n, n = 0, 1, 2, • • • . (70) 
Applying the condition (TTDl to (|67|) leads to the purely imaginary quasinormal mode as 



47rT> 



+ 



1 + A/1 + 



2 2 L 2 



(71) 
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which is similar to (I56"j) of the AdS2 black hole. In the large n limit, the QNFs behave as 



k 2 =0 



4nT, 



-in, 



11 



(72) 



which will be used to derive the area quantization of 4D2 Lifshitz black hole. 

Finally, to confirm the QNFs of s-mode, we start with equation ( 159)) with k 2 = 0, 

frn 3 + -, {y^ 2 -r%L 2 m 2 ){\-y^ _ 

2 y(i-y*) 2 + rMi-z/ 4 ) 2 2 

The solution to (1731) is given by the hypergeometric functions 



(73) 



H 2 (y) = d iy '- a (l-y*y* 2 F X 



2 



a 1 8 a 1 

- + -, - h -, 1 

4 2' 2 4 2' 



a 



+ d 2 y 2+a (l-y i yi 2 Fi 



B a 1 
— h u - 

2 + 4 + 2' 



a; 



/3 a 1 

2 + 4 + 2' 1+ 2 ; V 



(74) 



with c?i and d 2 arbitrary constants. At infinity (y — 0), c?i should be zero because of 2 — a < 0. 
Near the horizon (y —> 1), using the connection formula [35J, the solution ( I74p becomes 



H 2 (y) = d 2 y 2 + a Ci(l - yT> + C 2 (l - l/ 4 )^ 



(75) 



with the coefficients 



r(i + f)r(/3 



r(i + f)r(-/3 



(76) 



Near the horizon the first term in (1751) corresponds to the outgoing mode, while the second 
term corresponds to the ingoing mode. To have the ingoing mode near the horizon, we impose 
Ci = to give 



a 8 1 

- + - + - 
4 2 2 



-n, n = 0, 1, 2, • • • . 



The condition of ( 177)) provides us purely imaginary QNFs as 

,2 



-ir. 



4n + 2 + V4 + L 2 m 2 



(77) 



(7f 



which is exactly the same form as in (ITT)) . 

Consequently, we have shown that s-mode QNFs of 4D2 Lifshitz black hole are purely 
imaginary. Our question is " Can QNFs of Lifshitz black hole remain purely imaginary 
form even if k 2 7^ ?". At this stage, we could not answer to this question. However, the 
observation of QNFs of Lifshitz black holes suggests that /c 2 -term includes as a term in uj 
with ojr = 0. 
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5 Area spectrum of Lifshitz black holes 



It was suggested that black holes could provide a test bed for any proposed scheme for 
quantum theory of gravity. In this direction, Hod has combined the perturbation of black 
holes with the quantum mechanics and statistical physics to derive the quantum of the black 
hole area spectrum [29|. For a highly excited Schwarzschild black hole, Hod has used the 
real part ur of QNFs to obtain the area quanta of AA n = 41n[3]Zp. However, it is not 
consistent with AA n = 8nlp which was obtained from the fact that the black hole area is 
adiabatically invariant by Bekenstein [31J . Kunstatter has shown that the area spectrum is 
equally spaced for higher dimensional Schwarzschild black holes when using the adiabatically 
invariant integral [12] 

dE - (79) 



u(E) J uj r 

where (E, u>) are (energy, vibrational frequency) and (M, ur) are (black hole mass, real 
part of QNFs). On later, Maggiore has proposed that a black hole perturbed by external 
field is considered as a collection of damped harmonic oscillators [30J. He has regarded 
tu = ^Juo 2 R + uj\ as a physically proper frequency and thus, u = uj was used to derive 
AA n = 8irlp for highly excited QNFs of uj ^> lor by considering the transition from n to 
n — 1. 

Since QNFs of Lifshitz black holes are known to be purely imaginary and all thermody- 
namic quantities of Lifshitz black holes are known, we could use these to derive area spectrum 
of Lifshitz black holes by using Maggiore's method solely. In other words, the Hod's method 
is inapplicable to extracting information on Lifshitz black holes. Especially, we have their 
asymptotic frequencies which are given by (120"]) for 3D Lifshitz black hole, (153]) and (172|) for 
4D Lifshitz black holes. We review how to derive the quantization of horizon area for the 3D 
Lifshitz black hole [25]. We compute the adiabatic invariant I 

dM r+ 

(80) 



2G 3 (\/6-2)' 



where 



47rT H (V6-2) =2(V6-2) -± (81) 



Using the Bohr-Sommerfeld quantization condition of / ~ nh and considering the horizon 
area A = 2iir + , the quantized area spectrum is given by 

A n = AttG 3 (V6 - 2)nh, (82) 

which is not the universal area spectrum of A^ = 8irnh. The entropy spectrum takes the 
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form 

Sbh = ^= MVQ - 2)nh. (83) 

G3 

For 4D1 Lifshitz black hole, the adiabatic invariant / is given by 

f dM - V2T + (HA) 

-J— -4^G? (84) 



<4D1 



where 

U C = l [(wSl)n - (wSl)n-l] = 2ttTh = r 2 . (85) 

Considering the horizon area A = V 2 r\, the quantized area spectrum is given by 

A n = AnG 4 nh, (86) 

which is the universal area spectrum of A™ = 87rnh with G 4 = 2. In addition, the entropy 
spectrum can be obtained by 

5 4 " D1 = A = ^ (87) 
For other 4D2 Lifshitz black hole with L 2 = 1, the adiabatic invariant I takes the form 

dM V 2 r\ 



'4D2 



(88) 



where 



(^fcLo)n - (^£L ) n -l 



4ttT h = 4r 2 . (89) 
Considering the horizon area A = V 2 r+, the quantized area spectrum is given by 

A n = 8nG 4 nh, (90) 

which is not the universal area spectrum of A™ = 8imh with G 4 = 2. In this case one finds 
that the entropy spectrum is given by 

S2 D2 = A = 2nnh. (91) 

6 Discussions 

First of all, the purely imaginary QNFs show that scalar perturbation has no considerable 
oscillation stage around the Lifshitz black hole. This implies that the equilibrium is stable 
and thus, it is difficult to deviate the black hole from its equilibrium configuration. If the 
Lif/ CFT correspondence exists really, the thermalization timescale of the boundary conformal 
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field theory is given by r = A J T which implies that at high temperature, the field theory time 
scale is very small, indicating that a perturbation in the boundary conformal field theory is 
not long-lived and it decreases exponentially to zero. 

We regard the purely imaginary QNFs as an interesting feature of the Lifshitz black hole, 
in comparison to the AdS black hole. In this work, we have connected it to thermodynamic 
properties of Lifshitz black hole. All heat capacities of Lifshitz black holes are positive and all 
their free energies are negative, which means that all Lifshitz black holes belong to the GSBH. 
These globally stable Lifshitz black holes provide purely imaginary QNFs when choosing a 
scalar perturbation, which means that it is hard to take the black hole out off the equilibrium. 

Although we have tested a few of Lifshitz black holes, we suggest that most of Lifshitz 
black holes provide purely imaginary QNMs even if one uses different physical field (being 
not a minimally coupled scalar) as a perturbation field. 

As a byproduct, we have computed area spectrum of Lifshitz black holes by using the 
Maggiore's method. The constant 7 (or a) of A n = jnh (S n = anh) is given by 7 = a = 
4( v / 6 - 2)tt for 3D Lifshitz black hole with G 3 = 1, 7 = 8tt (a = vr) for 4D1 Lifshitz black 
hole, and 7 = 167r (a = 2tt) for 4D2 Lifshitz black hole. Even though all are not universal 
area spectrum, the area quantum spectrum could be derived from QNFs of Lifshitz black 
holes. Consequently, this result shows that both the horizon area and the entropy can be 
quantized for the globally stable Lifshitz black holes. 

We conclude with a comment on the recent issue related to scale covariant metric [43] 
which leads to a violation of hyperscaling of the dual field theory [H]. When a hyperscaling 
exponent 9 is zero, a scale covariant metric reduces to a scale invariant metric of Lifshitz 
metric. It would be interesting to answer to the question on the possible computation of 
QNFs of general Lifshitz black holes with scale covariant metric. Unfortunately, we could not 
obtain an analytic solution to the Klein-Gordon equation in these backgrounds (e.g., Eq.(5.2) 
of the X. Dong et al. work [43J). However, according to the AdS/CFT correspondence, it is 
known that the poles of the retarded Green's function in the momentum space correspond to 
QNFs [15]. If one computes such poles for the general Lifshitz black holes using the Lif/CFT 
correspondence, it might shed some light on finding QNFs. This is surely beyond the scope 
of the present paper, but nevertheless it is worthwhile to be explored in future work. 
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